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ABSTRACT 

The classical model of the chromosphere-corona transition does not account for the 
observed behavior of the differential emission measure for T < 10 s K. Several models have been 
proposed to resolve this discrepancy in physically different ways. Because the observed 
differential emission measure at the prominence-corona interface is on average nearly the same 
as in the chromosphere-corona transition, prominences offer a fresh testing ground for models 
tailored to the chromosphere-corona transition. 1 consider three such models and conclude 
that none extends in a natural way to the environment of prominences. I advance a simple idea 
involving thermal conduction both along and across the magnetic field from the corona into 
cool threads. 


1. Introduction 


The transition between chromospheric (or prominence) plasma and coronal plasma may be 
broadly defined as all gas in the temperature range 10 4 < T < 10 6 K. Throughout most of this 
range ( T > 3xl0 4 K), the plasma radiates mainly through optically thin permitted transitions in 
ions that exist with significant abundance only in relatively narrow intervals of temperature 
(AlogT < 0.3). For such a transition, the power radiated from a volume of plasma is 


P = ft J C(T) n e 2 dV « ft G(T max ) J n 2 dV . 


(1) 


Here C(T), which expresses the relative population of the upper level, is sharply peaked near 
7 'ma X and ft is an amalgam of atomic parameters and abundances. 

It is clear from this that the basic information observations provide about the thermal 
structure of the transition plasma is contained in the behavior of the differential emission 
measure (DEM), 


<(T) = n e 2 (T) 


din T ’ 


( 2 ) 


where dV(T) is the volume of plasma in the logarithmic temperature interval dlnT. 

In Figure 1, curve A represents the “average” behavior of the solar DEM, as compiled from 
many sources in the literature. The precise shape of the curve should not be taken too literally. 
There are variations from region to region, the spectroscopic diagnostics effectively smooth the 
curve to a resolution AlogT »0.3, and it is difficult to quantify all the sources of uncertainty in a 
given determination of the DEM. However, the basic shape is well established: a broad 
minimum in the range 5.0 < logT < 5.4 with a steeper rise to lower temperatures than to higher 
temperatures. The same curve can characterize the spatially-averaged emission from quiet 
regions, active regions (if the curve is shifted upward by about a factor of ten), and coronal 
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holes (cf. Chambe 1978; Raymond and Doyle 1981; Dere 1982). For the present purpose, the key 
observation is that this property extends to the spatially-averaged emission from prominences: 
the DEM curve of the prominence-corona interface is essentially similar, in shape and in 
strength, to the DEM curve of the quiet chromosphere-corona transition (Schmahl and Orrall 
1979; Yang, Nicholls and Morgan 1975). 1 Below, we explore some of the consequences of this 
uniformity. 


2. Models of the Chromosphere-Corona Transition 

In the standard picture (Giovanelli 1949; Athay 1966), the thermal structure of the transi- 
tion is determined mainly by thermal conduction along the magnetic field . The temperature 
dependence of classical thermal conductivity, k(T) oc T 5 / 2 , leads in a plane-parallel model to the 
familiar thin transition layer with a steep temperature gradient. The cooler part of the transi- 
tion (4.3 < logT < 5.0) is so thin (< 200 km) that it cannot radiate enough to account for the rise 
in DEM below T ~ 10 5 K. 

The standard model has been refined to allow for nonplanar magnetic geometry and steady 
flows. The result is that the basic thermal structure is little changed (Athay 1981, 1982). To 
date, no variant of the standard model has reproduced the observed low-temperature rise in the 
DEM. This failure has led to several recent models, three of which are next described. 

The three models share one feature: the cooler transition plasma occupies volumes 
separated in space or time from the volumes containing hotter transition plasma, and the hot 
branch of the DEM curve (logT > 5.3) is assumed to arise in the manner described by the stan- 
dard model. Otherwise the models are quite different. Because each (by design) can reproduce 
the cool branch of the DEM curve for the chromosphere-corona transition, other information 
must be brought to bear — e.g., observations of velocities, time dependence, or spatial structure 
(or other stars; see Section 4). Here we shall be concerned with the similarity between the DEM 
curves for the chromosphere-corona transition and the prominence-corona interface. Our 
brief treatment of the models is directed toward that end and is not intended to review all the 
advantages and disadvantages of each. 

a) Cool “ Coronal ” Loops 

The most recent model (Antiochos and Noci 1986) is in some ways the most natural. It 
also appears to be the least tenable as a model for the prominence-corona interface. 

In this model, the cooler transition plasma resides in low-lying magnetic loops that are 
nowhere hotter than a few times 10 s K (see Feldman 1984; Rabin and Moore 1984). A pleasing 
feature of the model is that cool loops arise as a solution of the standard “coronal loop equa- 
tions” on an equal footing with familiar hot-loop solutions. No additional mechanism is 
invoked to produce the cooler transition plasma. 

The equation of static energy balance, V-F c = H-C, where F c is the thermal conductive 
flux and H and C are the volumetric healing and cooling, admits an isothermal solution in the 
absence of gravity. To see this, one may picture the “loop” as a horizontal cylinder held at the 
same temperature on both ends (the chromosphere). The minor influence of gravity introduces 
a mild temperature gradient. Antiochos and Noci show that such near-isothermal solutions can- 
not exist if gravity is important, for then hydrostatic equilibrium dictates a rapid decrease in 
density with height. The radiative cooling, 

C(T) = n e 2 A(T) , (3) 

1 Systematic differences in detail do show up in a plot of the intensity ratio between prominence and quiet Sun for 
lines at various temperatures (Orrall and Schmahl 1976; Mariska. Doschek and Feldman 1979). 
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where A(T) is the radiative loss function, will also tend to decrease rapidly, whereas there is no 
reason to suppose that the heating is a similarly rapid function of height. Therefore, T must 
increase in order to increase A(T). However, beyond T « 10 5 K, A (T) becomes a decreasing 
function of T, and the demands of energy balance are incompatible with hydrostatic equili- 
brium. Therefore, the vertical extent of the cool loops must be less than the gravitational scale 
height at 10 5 K, about 5000 km. 

Antiochos and Noci show that the DEM curve produced by a single cool loop has negative 
slope, as is required by observations; but the value of the slope is sensitive to the characteris- 
tics of the heating function. The DEM from an assemblage of loops will also depend on the 
relative numbers of loops of different heights (different maximum temperatures) and the varia- 
tion of the maximum heating rate in a loop, F/ max , with the height of the loop. Antiochos and 
Noci use an ad hoc heating function parameterized as a function of magnetic field strength, 
which is in turn parameterized as a function of height. The heating function that gives a good 
match to the observed DEM is a steep function of loop height, H max oc /i loop -4 . In effect, the rise 
in DEM below 10 5 K is secured by preferentially heating the lowest loops. 

The sensitivity of this model to the distribution of assumed heating is unappealing, but the 
chromosphere-corona transition may nonetheless choose to meet the proper conditions. By the 
same token, however, it is hard to see why the heating in a prominence should conduct itself in 
the same way to produce the same DEM curve. 

There is a more fundamental objection to this model as applied to prominences. Quiescent 
prominences typically extend to heights of 30000 km or so, far beyond the limit of 5000 km set 
by the gravitational scale height at T ~ 10 s K. Yet emission from cooler transition plasma is 
observed over the entire prominence, not just near the footpoints. It is possible to imagine that 
cool-loop solutions exist high above the chromosphere by being trapped within vertical undula- 
tions of the magnetic field (Antiochos, personal communication). The undulations would have 
to have a characteristic peak-to-peak amplitude of less than 5000 km and would have to be an 
essential feature of prominence geometry. In the absence of independent evidence, this must 
be considered artificial. 


b) Heating-Cooling Cycles in Spicules 


This model (Athay 1984) builds on two observations. First, spicules usually disappear 
from view while still rising, presumably because they are being heating and get too hot to radi- 
ate H«. Second, a mean downflow is observed at T~10 5 K, and the mass flux is approximately 
the same as the upward mass flux in spicules (which must in any case return, since it is far 
greater than the mass flux in the solar wind). The model envisions three types of regions in the 
solar atmosphere: heating up; cooling down; and near equilibrium. It is assumed that the latter 
two regions are adequately described by conventional models of the chromosphere-corona 
transition and therefore suffice to produce the hot branch of the DEM curve. The cooler branch 
arises from spicules in the process of heating. 

The principal attraction of this model is that it explicitly recognizes the kinematic and 
time-dependent nature of the atmosphere. The principal disadvantage is that the slope of the 
cooler branch of the DEM curve is sensitive to the (unspecified) nature of the heating. 

In the simplest nonstationary model, with thermal conduction ignored, the DEM is deter- 
mined by how long a given volume of gas remains in each interval of temperature, which is in 
turn determined by the excess of heating over cooling, 


dT x T 
dt ~ 4.5 p 


( H-C ) , 


(4) 


where p is the gas pressure (taken to be constant during the heating). This' leads to the follow- 
ing expression for the DEM: 
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(5) 


c(T) = 9k(n e T) 2 F p j ~Zq 

where F p is the upward proton flux. 

It is evident that the emission measure directly reflects the temperature-dependence of the 
heating, <(T) oc H(T)T~\ except at temperatures for which H{T) only slightly exceeds C(T). 
Athay (1984, Fig. 1) secures agreement with observation by choosing H = constant and H = 
1.07 C at logT = 4.5 (the lowest temperature considered). Because C(T) (at constant pressure) is 
also nearly constant in the range 4.5 < log 7 < 4.9, the behavior of the small difference H-C is 
able to influence the DEM (steepening it, since C(T) decreases slightly). Note that both a close 
match between H and C at logT = 4.5 and the constancy of H are essential to this result. The 
match at logT = 4.5 itself presents a problem. Athay restricts his analysis to temperatures 
above the peak of La radiation at logT ~ 4.25 on the reasonable grounds that not all La radia- 
tion need arise in spicules and that radiative transfer in L« should not be neglected. However, 
unless one begs the question of how spicules are heated to logT = 4.5, the heating rate must be 
sufficient to surmount the La cooling peak, which is about 50 percent stronger than the cooling 
at logT = 4.5. If H is so chosen, H-C (for constant H ) varies by less than a factor of two over 
the range 4.5 < logT < 5.0. 

In summary, this model agrees with observation only for special choices of the heating 
rate. As in the model of Antiochos and Noci, this cannot be ruled out. But it is also clear how 
remarkable it would be if this attractive heating-cooling picture, inspired by spicules, were not 
only to carry over qualitatively to structures as different from spicules as are quiescent prom- 
inences, but were to be characterized by the same quantitative behavior of the volumetric heat- 
ing. 


c) Heating by Filamentary Electric Currents 


The key to producing the cooler branch of the DEM curve in this model (Rabin and Moore 
1984) is the temperature dependence of cross-field thermal conductivity. At constant pressure, 

K ± OC T~ 5 /2 (6) 

whereas 

CC T+ s /2 ; (7) 

so it is not surprising that cross-field conduction leads to a negative slope in the DEM. 

The model postulates that the cool transition plasma is internally heated by field-aligned 
electric currents. The current density must be high, 

j > 10 7 statamp cm -2 , (8) 

in order to balance radiative losses, but the current filaments (or sheets) must be thin, 

d < 1 km , (9) 


or else the magnetic field associated with the current itself will be unreasonably strong ( B > 
100 gauss). If the characteristic thickness of the filament is less than the cross-field conduc- 
tion length. 


d<d c = 


' «iT 
, n 2 A 


1/2 

~ 0.1 B _1 km , 


( 10 ) 


then cross-field conduction dominates the thermal structure, and the DEM curve has a slope in 
good agreement with observations [see equation (11) and Fig. 1]. 

None of the models considered here physically incorporates the heating mechanism. The 
previous two models simply include it as a parameter. The filamentary-current model specifies 
the mechanism (ohmic dissipation) but does not explain how the necessary currents are 
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produced. In the context of prominences, the present model has the further shortcoming that it 
does not predict the amount of cooler transition plasma. That is. although the shape of the 
cooler DEM curve is specified, its magnitude at a given temperature, and therefore the tempera- 
ture at which it crosses the hot branch of the DEM curve, is determined by the number of 
current filaments per unit area. Why should this number be comparable in low-lying loops and 
in quiescent prominences? [However, as illustrated in Fig. 1, the steepness of the cool branch 
of the DEM means that it may be substantially displaced in magnitude (vertically) without much 
changing the temperature of its intersection with the hot branch.] 


3. Thermal Conduction Across and Along the Magnetic Field 


It is instructive to analyze the models considered above because each has interesting (but 
different) physics and each has a tantalizing set of strong and weak points. Yet, in the final 
analysis, the close correspondence between the DEM curves for the chromosphere-corona tran- 
sition (which itself encompasses everything from active regions to coronal holes) and the 
prominence-corona interface serves as much as anything to draw attention to the most attrac- 
tive feature of the original conception of the transition region: that thermal conduction “does 
everything” — provides the energy, determines the thermal structure, and is there whenever 
and wherever needed. In contrast, the cool-coronal-loop model and the healing-spicule model 
rely on the detailed behavior of an unspecified heating mechanism to match the cool branch of 
the DEM — they must be “fine-tuned.” The filamentary-current model requires that the current 
elements be present, in comparable numbers, in widely different environments. 

One possibility for recovering a unified treatment of the transition plasma is that nonclas- 
sical electron transport influences the structure of the transition (Roussel-Dupre 1980a, b\ 
Shoub 1983). It appears that streaming by fast electrons in the enhanced (non-Maxwellian) tail 
of the distribution function has only a minor effect on the excitation and energy balance of ions 
used as transition-plasma diagnostics (Keenan 1984; Owocki and Canfield 1986). The impor- 
tance of the nonclassical heat flux, both in terms of magnitude and distribution over tempera- 
ture, is not yet clear. The nonclassical flux depends sensitively on the temperature structure of 
the transition, which has thusfar been assumed rather than calculated self-consistently (Owocki 
and Canfield 1986). 

Here 1 consider a different idea, not yet developed to the stage of a “model.” Any cool 
magnetic structure that protrudes into the corona will receive heat both along and across field 
lines. 2 Classical cross-field conduction is usually ignored at coronal temperatures because it is 
much weaker than along-field conduction. At T ~ 10 5 K, the disparity is smaller but still sub- 
stantial. An analytic model with planar geometry, including conduction and radiation (in 
power-law segments) gives the following expressions for the DEM per unit area exposed to the 
corona: 


<±(T) = 10 227 (n e T) l5 2 B- 1 T, ow4 3 / 2 [l-(T low /T) 3 
<ii CO = 10 262 (n e T) ls [l-(T low /T) 1 /2j~ 1/2 T 5 5 /< C m- 5 , 


Tc~ 7 ! 2 cm' 


(ID 

(12) 


where T low is the cool temperature at which the conductive flux goes to zero and T 5 means T in 
units of 10 s K. The magnitude of tn(T) is about right to match observations at T = 10 5 K. 
Therefore, to turn up the DEM below 10 s K, we must increase the area exposed to cross-field 
conduction. In Figure 1, curves B-D show the results of simple area— weighted superpositions 
with the cross-field area increased by factors of 10 35 , 10 40 , or 10 45 (the curves have been dis- 
placed upward one decade from the mean observed curve). Clearly, a typical ratio ~10 4 is 


7 1 

Appropriately, it seems to have been in the context of prominences that cross-field thermal conduction was first con- 
sidered in the solar atmosphere (Orrall and Zirker 1961; Doherty and Menzel 1965). 
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DIFFERENTIAL EMISSION MEASURE 


Fig. 1. Observed anc! computed DEM curves. 
A: “average” observed curve for non-active 
areas; a mean curve for active regions would 
be higher by about a factor of ten. B-D: 
area-weighted superpositions of transverse 
and longitudinal thermal conduction, with 
the area exposed to cross-field conduction 
increased by factors of 10 35 , 1 0 4 °, and 10 45 
respectively; the curves have been displaced 
upward by one decade. E-G : mixtures of 
transverse and longitudinal conduction for 
magnetic field inclined to the isotherms by 
0.05, 0.1 and 0.2 degrees respectively; the 
curves have been displaced upward by two 
decades. In the computations, n e T = 10 15 , 
r low = 10 4 , and B = 3. 


4.0 4.S 5.0 5.5 6.0 6.5 

LOG TEMPERATURE 



needed to match the observed curve, but, because the cool branch of the DEM is steep, the ratio 
can vary by more than an order of magnitude from place to place without moving the minimum 
of the DEM curve out of its observed range, log 7" ~ 5. 0-5.4. It may be noted that the slope of 
the hotter branch of the DEM in curves B-D (as well as E-G ) is shallower than the observed 
slope. This is a consequence of the simplified radiation law and planar geometry adopted in 
the present treatment, which does not attempt to model the hotter transition per se {cf. Athay 
1981, 1982). 

If magnetic field lines characteristically leave cool structures at a small angle to the iso- 
therms, transverse conduction and longitudinal conduction are physically blended rather than 
superposed. By rotating the conductivity tensor, one may still solve the planar problem analyti- 
cally. The result is shown for three angles by curves E-G in Figure 1; the curves have been dis- 
placed upward by two decades. The shape of the observed curve is reasonably well-matched 
for an angle of 0.1 or 0.2 degrees. As in the previous example, the area exposed to the corona 
must be increased by about four orders of magnitude in order to reach the proper normaliza- 
tion. 

The necessary area exposed to cross-field conduction will be smaller if the cross-field 
transport is nonclassical. For Bohm-type diffusion, the emission measure is increased over its 
classical value by 


— * 10 08 B>/2 (n e T) 15 ~ 1 / 2 T 1ow4 -s/ 4 T 5 s /2 , (13) 

<dc 

about an order of magnitude at T = 10 5 K. However, the slope of the DEM curve produced by 
Bohm conduction is approximately -1, shallower than the observed curve. And, where the 
wave-particle interactions that lead to Bohm-like conduction are active, parallel thermal con- 
duction will be inhibited (and changed in temperature dependence) at the same time that per- 
pendicular conduction is enhanced. 
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4. Discussion 


Although they are too simple, the analytic models considered in the previous section illus- 
trate the basic idea: thermal conduction can determine the slopes of both branches of the DEM, 
without adjustable parameters. The principal attraction of a single microphysical mechanism Is 
that it can operate basically unchanged in a variety of settings, from prominences to plages. 
Other stars may now be added to the list of environments. Observations with 1UE and with the 
Einstein Observatory have produced DEM curves broadly similar in shape (although not neces- 
sarily in strength) to the solar curve among stars from a range of spectral types (late F to early 
M) and gravities ( e.g ., Zolcinski et al 1982; Linsky et al. 1982). This again suggests a process 
without too many “knobs.” 

The price exacted for the idea considered here is a willingness to contemplate transition 
plasmas that are highly fragmented, such that the area exposed to the corona is three or four 
orders of magnitude greater than the projected surface area. Is this reasonable? 

Historically, whenever the resolution of our observations has improved, we have seen finer 
structure (prominences are a good case in point). There is no observational evidence to suggest 
that this process is nearing an end. There is far to go from the internal viewpoint of the transi- 
tion plasma: its fundamental scales (such as Debye length and proton gyroradius) are often cen- 
timeters or less. We must rely on indirect diagnostics. 

Schmahl and Orrall (1979) have presented evidence that Lyman continuum absorption 
affects EUV spectra everywhere on the Sun, from network cells to active regions to quiescent 
prominences (also see Withbroe 1977; Kanno 1979; Nishikawa 1983). The inferred attenuation, 
about a factor of six, is remarkably constant from place to place. Orrall and Schmahl (1980) 
have analyzed the Lyman continuum data in more detail for nine hedgerow prominences. They 
find that there must be at least 4-10 cool threads or sheets along a typical line of sight, but an 
upper limit is not determined. Further information on the distribution of absorbing material 
may come from the observation that the Lyman continuum absorption is smaller at 10 6 K than 
at 10 s K (Schmahl and Orrall 1986). 

Recently, Fontenla and Rovira (1985) have constructed NLTE models of prominence threads 
and have computed profiles and absolute intensities of La, L 0 , Ha and the Lyman continuum 
from an ensemble of threads (observations are discussed by Vial 1982). They infer that the 
minimum number of threads along a line of sight ranges from 10 to over 100, depending on the 
model; but again, no upper limit is placed on the number of threads. 

Such studies show promise for constraining the degree of filamentation in prominences. 
Already it is clear that imagining thousand fold filamentation along the line of sight is not as 
far-fetched as it might first appear. 

There are good opportunities for progress in at least two other directions. First, we need 
spatially and temporally resolved (~ 1", ~ 1 min) determinations of the DEM in various struc- 
tures. It will not be easy to collect enough photons! Still, it is vital to progress beyond “aver- 
age” DEM curves now that several models vie to reproduce them. The fact that Skylab and SMM 
observations sometimes show differences in spatial structure between various transitional and 
coronal lines should alert us to expect dramatic departures from the mean curves (Cheng 1980; 
Poland and Tandberg-Hanssen 1983). It may be possible to test whether the DEM depends on 
pressure and magnetic field strength in the manner indicated by Equations 1 1 and 12. 

A second fruitful direction will be further analysis of stellar DEM curves. Here the com- 
plete spatial averaging is offset by the opportunity to sample a wide range of coronal and chro- 
mospheric parameters. For example, n Dra (G2 lb-11), which shows strong chromospheric and 
transition emission in conjunction with a relatively weak and extended corona, poses difficulties 
for models of the transition plasma based on thermal conduction and points toward internal 
heating (Brown et al. 1984). 
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